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MATHEMATICAL  ASPECTS  OF  ELECTROMAGNETIC  THEORY  III 


I,  Bessel  Functions 

In  this  section  we  shall  investigate  the  Bessel  functions 
by  studying  the  solutions  of  the  two  dimensional  wave  equation 
written  in  polar  form.  N o  prior  knowledge  of  Bessel  functions 
will  be  assumed;  their  properties  ana  representations  will 
arise  as  a  natural  consequence  of  the  by  now  familiar  methods 
of  solution. 

The  wave  equation  in  rectangular  coordinates, 

^xx  +  ‘Pyy  +  ^  =  "  ^{x  "  xo}  <f(y  "  yo}' 
is  simply  a  generalisation  of  equation  (33)  Technical  Report 

No.  2  and  its  solution  is  easily  seen  to  be  equation  (34) 

written  in  a  slightly  more  general  form: 


.  .  1  {  ‘P(x-*0>  +  Wit2  -  p2|y-y  | 

^  m-&  W - 1— — - —  dp  , 

C 


where  the  contour  C  is  identical  with  the  contour  used  in  (34). 

However,  the  wave  equation  may  be  solved  in  any  coordi¬ 
nate  system  which  permits  the  separation  of  variables.  Thus 
we  can  introduce  plane  polar  coordinates  and  solve 


<prr+  r^+  T^ee  +  ***  * 


-  Ar-ro)  f(  9-0o) 


2 


In  order  to  accomplish  this  we  will  try  a  solution  of  the 

oo 

form  <p(r,0 )  =  un(r)ein9.  This  expression  for  <p(r,0) 

-co 

implies  that,  for  m  =  0,  1,  2, 

■  C  &  ’(r'e)»-lne  - 
*0 

also,  integration  by  parts  shows  that 

^  <P60e  d0  =  (<pQ  +  in<p)e“ln®  ^  n2q>e'in0d0  *  -n^(r)  , 

since,  from  physical  considerations,  we  can  assume  <p(r,0)  - 
=  ip(r,2it)  and  cp  (r,0)  =  <p  (r,2it).  Thus,  when  the  wave  equa- 
tion  is  multiplied  by  e"in0  and  both  sides  are  integrated 
over  0  <  0  <  2n,  we  have 


(1)  un  +  7uA  +  (k2  "  4,un  s  "  fii.  e'1 


,  0  <  r  <  oo , 


There  is  evidently  only  one  boundary  condition  for  equa¬ 
tion  (1),  namely,  un(r)  must  satisfy  the  radiation  condition 
and  be  outgoing  at  infinity,  A  boundary  condition  at  r  *  0 
is  necessary  in  order  to  define  a  unique  olution.  Mathe¬ 
matically,  the  change  to  polar  coordinates  causes  the  origin 
r  *  0  to  become  a  singular  point  of  the  differential  equation  (1 
As  a  consequence,  not  all  the  solutions  of  (1)  are  regular 


at  r  =  0.  But  physically  the  origin  is  a  point  in  space 


without  special  distinction,  and  the  field  produced  by  the 
line  source  along  (r  ,fl  )  should  be  defined  here.  There- 
fore  our  second  conditi  n  is  that  the  solution  of  (1)  which 
we  select  must  have  no  singularity  at  the  origin. 

The  usual  Green's  function  technique  will  be  applied  to 
equation  (1),  The  homogeneous  equation  is 

w"(r)  +  -pw^(r)  +  (k2  -  %)wn(r)  *  0, 

r 

By  setting  kr  =*  0 ,  we  obtain  Bessel's  equation 

(2)  wjip  )  +  ±w#p)  ♦  (1  -p)w„(p)  -  0, 
where  the  primes  now  denote  differentiation  with  respect 
to  p  ,  Attempting  a  power  series  solution  of  the  form 
w  *  f>a(l  +  ajf  +  a2p2  + 

we  find  after  substitution  into  (2)  that  the  indiclal  equation 
2  2 

is  a  -  n  *0,  i.e,,  a  s  +  n,  Thus  the  two  fundamental  solu¬ 
tions  of  equation  (2)  behave  as  pn  and  ^>“n,  The  solution 
that  behaves  as  will  be  discarded  because  it  is  not 

regular  at  p  =  0;  the  other  solution  will  be  labeled  Jn 
n  >  0,  and  called  the  Bessel  function  of  order  n.  Conse¬ 
quently,  for  r  <  rQ,  ufl(kr)  «  wfl(kr)  =  const  Jfl(kr )  <v/(kr) I 
denotes  "behaves  as"). 

Notice  that  for  n  *  0  it  appears  the  homogeneous  equation 


has  only  one  solution;  however,  in  this  case  It  can  be  shown 
that  a  second  solution  exists  which  behaves  like  log  p  • 

In  any  case,  the  equation  (2)  always  has  two  linearly  inde¬ 
pendent  solutions. 

In  order  to  determine  which  solution  of  equation  (2)  is 
outgoing  at  infinity,  we  make  the  substitution  wn  =  p~  '  vn» 
this  will  remove  the  first  derivative  term  in  (2),  (Such 
a  removal  can  always  be  accomplished  in  a  second  order  diffe¬ 
rential  equation  by  means  of  the  general  substitution  w  *  fv, 
where  f  is  a  suitably  chosen  function.)  The  homogeneous  equa¬ 
tion  becomes 

+  (I  -  - 0 

Since  we  are  interested  in  vn(p)  for  large  values  of  f>  , 

2  A  +  IP 

we  neglect  the  term  - — 5  and  thus  obtain  v  e” 

-1/2  +  1 if 

and  e  for  p  large.  The  outgoing  and  incoming 

solutions  of  equation  (2)  are  called  the  Hankel  functions 
and  defined  as  follows:  H^1  ^ (/3) 

where  n  >  0  •  Recall  that  the  factor  ~  has 

also  appeared  in  earlier  investigations  of  the  behavior  of 
cylindrical  waves  at  large  distances  from  their  source  and 
that  this  factor  is  compatible  with  the  conservation  of  field 


energy  (e,g.,  see  section  VIII  in  Technical  Report  No,  2) 


Proceeding  with  the  Green's  function  technique,  we  put 
ufi(r)  *  const  Jn(r),  t  <  rQ 

=  const  H^(kr),  r  >  rQ  . 

When  we  cross  multiply  in  order  to  assure  the  continuity 
of  un(r),  we  find 

Vr)  ■  enHi1,‘kro>Jn(kr)>  r<ro 

*  enJn(llrcX1>(kr)>  r>ro- 

The  constants  c_  are  to  be  determined  by  the  jump  condition, 

”  -ineo 

a.(r)|r  -u.(r)|p  - 

o  o" 

i.e..  the  c  should  be  such  that 

-Ine0 

(3)  Jn(kr)]r.r  -  . 

o  o 

If  we  let  the  prime  denote  differentiation  with  respect  to 


kr,  equation  (3)  can  be  rewritten 


c  k[J  ] 

n  n  n  n  n  r=r 


The  bracketed  expression  is  essentially  the  Wronskian  of 


the  two  chosen  solutions  of  the  homogeneous  equation.  It 

should  not  be  surprising  thatfrJ-H^1  ^  -  rH^^J']  * 

r  n  n  n  n  r*r 

-in©  ° 

*  const  e  Is  independent  of  rQ,  for  if  equation  (2) 

is  written  in  self-adjoint  form. 


we  find 


-  6  - 

2 

(p  "A>’  +  (/>  -pwr,  -  o. 

(/’JJ)'  -  If  -,-2>Jn-c’ 

and 

(  +  (p  -  jtn < 1  >  =  o 

=*  -Vpii11')'  -0 

=*  Vp^1’-/3  h"1|,j"1  =0 

^  P  IJn  HiU  -  Hn1)'Jn1  "  const- 

Because  the  Wronskian  is  a  constant,  the  next  step 

would  ordinarily  be  to  evaluate  Jn,  ,  and  their 

derivatives  at  some  convenient  value  of  r  ,  substitute 

o 

back  into  equation  (3),  and  thus  find  the  cn.  Usually  rQ 

is  taken  large  since  the  asymptotic  behaviors  of  the  Bessel 

and  Hankel  functions  are  well  known  and  easy  to  work  with. 

But  because  we  assume  no  knowledge  of  these  functions,  their 

behavior,  asymptotic  or  otherwise,  must  first  be  determined 

as  part  of  the  solution  of  equation  (1), 

To  begin  with,  recurrence  relations  can  be  obtained 

for  the  Solutions  wn(kr)  of  the  homogeneous  equation.  We 

know  w  (kr)ein®  satisfies  u>  +  <p,„,  +  =  0  :  we  notice 

n  txx  yy 

also  that  «  ,  <p  ,  and  any  linear  combinations  thereof  are 
x  y 

solutions.  Therefore, the  combinations 


+  1^>wn(Kr)®lr'e  <u,d  &  '  wn<kr,elne 


-  7  - 


satisfy  <p  +  <p  +  kfc<p  *  0.  If  and  are  expressed 

xx  yy  0  a 


-2. 


in  polar  coordinates  we  find  that 
)“»</,)  e‘ln*‘)9  and 

must  each  be  solutions  of  the  wave  eqviation.  Consider  the 
first  of  these,  which  reduces  to  (w^  - —WjjJe1  ^n+1  : 

by  referring  back  to  the  derivation  of  equation  (1),  it  is 
not  difficult  to  see  that  (w^  ”^wn)  solves  equation  (2) 
when  n  has  been  replaced  by  n+1  in  that  equation.  We  con¬ 
clude  that  for  some  constant  a, 

n  p  n  n+1  * 

Similarly,  using  the  second  of  the  linear  combinations 

above,  we  conclude  that 

w’  +  £w_  *  (3w_  , . 
n  p  n  n— l 

The  quantities  a  and  0  are  yet  to  be  specified;  to  this 

end,  note  that  w*  ,  +  ,  =  pw„,  from  which  follows 

n+l  p  n  n' 

u"  +  JL  w'  —  —jt  w  —  aOw  . 
n  pn  pZ  n  ^n • 

Because  wn  satisfies  the  homogeneous  equation  (2),  we  have 
apwn  =  and  aj3  =  -1.  We  shall  choose  a  *  -1,  (J  *  1* 

The  recurrence  relations  are  therefore 


(4) 


w »  -  ^w 
n  p  n 


-  wn+l  ,  Kn  +  jSwn  "  wn-l  • 


-  8  *• 

Recall  that  J  ,  and  are  solutions  of  the  homo- 

n'  n  '  n 

geneous  equation  and  so  are  qualified  for  the  above  relations; 
thus,  for  example,  if  can  be  determined,  the  above 

formulae  may  be  used  to  write  for  any  n.  This  fact  will 

enable  us  to  arrive  at  an  integral  representation  for 
as  soon  as  wefve  found  one  for  from  this  integral 

representation  will  follow  the  asymptotic  behaviors  of 
H^,  and  J  ,  which  in  turn  will  enable  us  to  compute 

the  cn. 

Since  Jn(kr)  (kr)^nl,  Jfl(0)  “  0  for  n  ^  0.  We  shall 

specify  that  J  (0 )  =  1  •  It  is  known  from  Section  VIII  that 

for  x=  y=  0,  the  solution  of 
o  o 

»xx  +  v  + k2‘f  “  -  £u-*o>  ^(y-y0> 

can  be  written 


(p  = 


const  \  eIltr  005  9de  . 


^  $  —  f>Car>4. 


y 


9 


On  the  other  hand.  If  r  3  0  we  have  seen  that 

o 


oo 

<P  =  TT  ^  oTlJn[WrJHi1\Kr)zlnB  n  >  r. 
n*  -oo 


o  'n 


Hence,  because  Jn(0)  =  0  for  n  /  0  and  JQ( 0)  3  1,  we  must  have 


H^lj(kr)  »  K  J  eikr  008  ede, 
C1 


where  K  is  some  constant.  Now,  H^^(kr)  r>s  (kr )“l/2eikr 

(kr)“1/2e“lkr,  which  implies  that 


and  H^2^(kr) 


H^2)(kr)  3  H*l,(kr)  3  K  J  eikr  coS  0d9  , 

We  wish  to  compute  the  complex  conjugate  of  the  integral, 
a  task  which  is  complicated  by  the  fact  that  lies  in  the 
complex  plane,  may  be  reduced  to  the  sum  of  the  three 

following  integrals  over  real  paths  by  the  substitution 
0  =  a  +  ip  : 


H*11'  -Hi  \  .tkr  c°»h  hf  +  K  j 
00  *0 


TOD 

eikr  cos  ada  +lk  (  #-lkr  cosh 


1 


therefore 


Mm 


hJ2L  inj  e"lkr  C08h  Pdp  +  k|  e~lkr  C08  ado  -  iKj  elkr  C08h  % 


11 


But  because  the  integrand  is  periodic  with  respect  to  a, 
the  integrals  taken  over  the  vertical  portions  of  cancel 
one  another.  Therefore, 


H^'lkr)  +  H*2)(kr)  = 


<2” 


ikr  cos  a. 
e  da  . 


The  above  sum  is  a  solution  of  equation  (2);  it  is  real 
since  ^  and  H^2^  are  complex  conjugates,  and  it  is 
regular  at  the  origin.  We  conclude  that  +  H^2^  is 

euqal  to  some  multiple  of  JQ;  we  define 

H^(kr)  +  H(2^(kr)  =  2J  (kr). 


From  the  fact  that  J  (0)  =  1,  we  find  that  K  =  it-*, 

o  ' 


We 


now  have  the  desired  integral  representations: 


(5)  H<1}(kr) 


■ii 


eikr  cos  9d9>  H‘2)(kr)  =  i  J  e 

c 


j. 


ikr  cos  © 


de. 


and 


jo(ki-) 


■if. 


0ikr  cos  aA„ 
e  da  , 


Note:  Since  H^.  h^,  and  J  all  satisfy  the  recursion 
■  ■  n  '  n  '  n 

relations  (4).  we  see  that  +  H^2^  =  2 J  ^  2J  . 

'  o  o  onn  n 

The  difference  of  the  two  nth-order  Hankel  functions  is 

usually  denoted  by  2N  =  ^  -  H^2^  and  is  known  as  the 

n  n  n 

Neumann  function  of  order  n. 

It  shall  now  be  shown  that  the  outgoing  Hankel  function 


of  order  n  is  given  by 


~  12  - 


(6)  h  yh/>) 


(-I)n  \  elf-  cos  eeinede  # 


'1 


Equation  (6)  will  follow  directly  from  the  fact  that  the 
integral  (5)  satisfies  the  recursion  formula  (!].)  for  Hn  , 
namely 

-h<‘>  =  h(‘)’  . 
n+i  n  ft  n 

In  order  to  prove  this  fact,  we  shift  the  path  of  integration 


eikr  cos  edQ 


eikr  COS  9de, 


"1 

and  a  minor  rearrangement  of  that  proof  also  gives 


elkr  cos  9eined0  =  ^eikr  cos  9eln9de  < 


Let  i=i£\e*  /,cos  Vned9  =  I£>n(f) 


(  p  =  kr);  then 


is  uniformly  convergent  for  §  <  £  <  it. 


since 


-  13  - 


I»  -  —i  =  {ri)H  ]  e1^008  ®*ln®[i 
£  ,n  p  £»n  it  J  9  l 

cfc 

,  t-‘.ln-  (etf  cos  Ggine 


cos  0  -  £]de 


z° 


(-i)n  e1/*008  ®e 

+  n 


in© 


[i  cos  e  -  sin  e]de  + 


where  we  have  integrated  by  parts. 


Since  the  first  term  here  equals  -1^  n+j  ,  it  need  only  he 
shown  that  the  integrated  part  is  aero.  But 


(i p  cos  e  ine 


+  e*  P  cos  vna 


s  eihcos(- ^  +  ip)ein(-^+  i0|°  + 

*  1 00 


+  #  +  ei/f  cos(it~E  +  ip)ein(n-f+ip)j 

-t  *"  * 


-  00 
0 


=  lim  [exp(~i  a  cos  ^  cosh  p  -  Z3  sin  €  sinh  0  +  in(n-t)  +  np)+,. 

p  -*  00 

...  +  exp{  i  p  cos(£-  ip)  -  ln£-  np)J. 

This  limit  is  equal  to  zero  so  long  as  it  >  £  >  0,  for  In 
the  first  term  sin£  sinh  p  +  np^  -  pslng  e^  +  np  —  -oo 
as  p  ■?  oo,  and  in  the  second  term  i^>  cos(fc  -  ip)  -  np  * 

=  i  f>  cos£  coshp  -f>  s  Inf  sinh  p  -  np  rv  - p  sin k  e^3  -  np  —  -cr» 
as  p -*  oo  .  Equation  (6)  Is  thereby  proved. 

Equation  (6)  enables  us,  after  some  familiar  prelimi¬ 
nary  manipulations,  to  work  out  an  asymptotic  expression 
for  If  we  let  p  *  cos  0,  (6)  becomes 


HJ'V 


-  14  - 

i=i£  (  eJ(0  Pelne 

v.  J  rr 


s/W?  ' 


p-plane 


i  + 6# 


- 1  ^  C 

As  before,  in  order  to  obtain  an  integrand  which  decreases 
as  p  increases,  the  contour  C  is  shifted  to  CQ  and  p  put 
equal  to  1  +  is  (notice  that  ein0  is  analytic  except  at 
p  =  +  1),  It  follows  that 

GO 

H<>V)  -  |«I<|5  ’  tA  ■  B‘/2)  (  jlfaiy  a.  . 

"  •  11  J  y/7  l/2^Is 

0 

Since  ein0  has  no  derivative  at  s  **  0  ein®  * 


=  inein0[  - 


=  me  _i  -  (Rspysaecws]  )  we  may  not  legitimately  use  the 
\J»  -  2is 

theorem  of  Section  VIII.  However,  at  the  end  of  that  section 
an  alternate  method  was  used  to  arrive  at  an  asymptotic 
approximation  of  E  ,  and  this  method  will  also  find  appli- 
cation  here. 

The  value  of  eln0(2  +  is)“  ^ 2  at  s  =  0  is  1»(2)~  1/2  • 


we  can  write 


1  ,  r  { 2  elT>9  -  <[2»ls  . 

l  I— «- 


V  2+is  {T  \JT\l2  +  is 

Substituting  this  into  the  expression  for  and  integrating 

the  leading  term,  we  have 


where 


./r.*v  -  ^  -  n*/2)  ♦  e. 


00 


We  may  write 

jr«tB9  -  m  .ln9  - 1  +  a  JST, 

{S  iSZ 7T  v  2+is  \|T  V2+is 

The  second  term  on  the  right,  which  also  appeared  in  Section  VII 
is  absolutely  bounded  by  — ,  A  bound  for  the  first  term 
can  also  be  obtained.  Recalling  that  cose  *  1  +  is,  it  is 

easy  to  see  that  (e10  -  1)  »  j2lsel0//2;  and  because 
(ein0  -  1)  *  (ei0  -  +  el(n-2)e  +  +  e10  +  j), 

it  will  follow  that  |eln0  -  lj  <  ^2  n  tfs^  Hence, 


This  error  estimate  is  satisfactory  whenever  n  « ^  (the 
behavior  of  for  large  n  will  be  investigated  later). 

Thus,  since  H<2)(^}  *  and  2Jn(f>)  -  H*1^)  +  H<2)(p), 

we  have  determined  the  asymptotic  behavior  of  the  Bessel 


and  Hankel  functions: 


16  - 


(Hn1>‘ "A'WZ) 

L(2)(o  y.-j  2  ( P  m  *A  “  n*/2) 

(7.uHn  ’ 

jn(p)^^^c°s(p  -*A  -  n«/2)  (f  >>  n  >  0) 

We  are  now  finally  in  a  position  to  compute  the  cn  of 
equation  (3)  and  thus  complete  our  solution.  In  general, 
asymptotic  forms  may  not  be  differentiated  to  yield  the 
asymptotic  forms  of  the  derivatives.  However,  in  this  case, 
since  the  paths  Cj  and  which  appear  in  the  integral 
expressions  of  and  H^2^  may  be  slightly  altered  in 

order  to  assure  the  absolute  convergence  of  the  integrals, 
equation  (6)  can  be  differentiated  with  respect  toy0.  The 
estimation  of  may  be  then  carried'  out,  and  by 

referring  to  the  derivation  of  the  asymptotic  form  of 
H n^(p)  4 4  is  almost  immediate  that 

~  +  OCp”1)  • 

Therefore, 

*A  "  I”'/,2)  *  IHn* 

w>  sin(p  -  n/l i  -  n«/£), 

where  the  terms  containing  have  been  dropped.  When 

these  asymptotic  forms  are  substituted  into  equation  (3), 
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j i  “ in0_ 

we  find  c  [  ~]  *  -e  . 

n  oo 

The  solution  <p  =  7  u  (r)ein®  of  the  wave  equation  is 
T  i-  ■  ”  n 
-oo 

00 

.  / « \  in( 0-0ft ) 

,flr'e)  "  k2l  HM(kro>  J|n|(kr)e  '  r<ro 

-00 

oo 

4  { i  \  in(e-e  ) 

*  Jtn|  (kro,HInf  (kr)e  *  r>  ro 

-00 

It  may  be  more  simply  written  by  defining 
r>  §  max(r,ro)  ,  r^  5  min(r,rQ): 
we  get 

jjSL.  /jj  in (9-8  ) 

(8)  *(r,e)  >_  Hjiflkr,)  J)n|(kr.>  °  . 

-CO 

*  *  * 


When  the  integral  expression  for  H^Ukr)  was  being 

c  (  2 ) 

computed  we  saw  that  for  rQ  =  0,  <p(r,e)  3  '(kr); 

_i  -ineQ 

c  is  known  from  c .  =  4r  «  •  Therefore,  if  (r  .0  ) 

o  n  c  o  o 

on  the  xy-plane  is  taken  as  the  origin  we  can  write 


16 


Equation  (9)  leads  to  the  evaluation  of  a  definite  Integral, 
Since  R  ■  \j  ro  +  f2  "  2rQr  cos(e-eo)  ,  it  follows  that 

2it 


cn  r - 

Ho1)(k  V'!  +  r2  -  2r0r 


co5(e-eJe“in8de  ® 


=  H(i) 


sHM<kr>>J|n|<kV« 


-ine. 


We  might  ask  what  happens  when  the  line  source  is  at 
infinity,  i.e,,  when  rQ  -*  oo  •  As  rQ  becomes  large, 

Rss  roC1+J7"  ~  ¥■  co»(®-e  n 1/2  = 

rQ  ro 

m  r°  \j  "  ^  cos{0-eo)]  +  0(  r"1 )  , 

and  referring  to  the  asymptotic  form  of  H^l)(kR),  we  have 

H‘uo*)  -<4r>1/2«1(kR  '  %A)  ♦  oi*'1) 

(U)  *  (  >§?  )1/5e1(ltr<J  -  *A).-*kr  co*(e-80>  +  0(r-l. 

*o  o  / 

-ikr  cos(e-e  ) 

Since  the  function  e  °  represents  an  incoming 

plane  wave  from  the  direction  0o,  we  see  that  the  line  source 
at  very  large  distance  from  the  origin  is  asymptotically 
a  plane  wave  of  very  small  amplitude  (r~2/2). 
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Consider  the  right-hand  side  of  equation  (9)»  Since 

r  is  fixed  and  r  —  cd  .  we  have  r  =  r  and  r^  *  r. 

o'  >  o  < 

With  the  aid  of  estimates  of  and  for  large  values 

n  n 

of  n  that  we  shall  obtain  later ,  it  is  easy  to  'Show  that 
the  series  in  (9)  is  absolutely  and  uniformly  convergent. 
In  fact,  we  can  show  that 


oo  N 


U2)  4  Hf^(kro)J|n|<1'r>‘lnl9'e°)+:X'-^'1> 

If  rQ  is  large  compared  to  N,  equation  (7)  gives 

»i:k>  •  (&.)-*  «-;■> 

Substituting  (11)  and  (13)  in  (9),  we  obtain 

00 


-ikr  cos(e-eJ 
(111-)  c  0 


in ( 0-0  ) 

(-i)nJ|n| (kr)e  0 


-oo 


a  representation  of  a  plane  wave  as  a  sum  of  incoming 
and  outgoing  cylindrical  waves. 


It  still  retains  to  discuss  the  asymptotic  behavior 

/ ,  >  /  i  \n  C  1  pcos  0  in© 

of  J  e  e  d@  for  the  case  of  large 

u  n  *  X 
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(the  path  is  shown  on  page  12)  •  Instead  of  however, 
we  shall  use  the  contour  C* : 


Note  that  -*/2  <  a  <  n/ 2  on  C',  a  choice  of  a  the  reason 
for  which  will  become  apparent  below.  Since  our  concern 
here  is  with  large  n,  we  ought  to  Include  n  into  the  change 
of  variable;  thus  in  the  integral  for  put 

p  3  f*  cos  e  +  ne.  This  results  in 


where  &  is  computed  from 

p  3  p^  +  ip£  *  cos  a  cosh  (J  +  na  +  i(-sin  a  sinh  0  +  np): 
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The  loop,  whose  length  depends  on  n,  arises  -  from  the  fact 
that  when  a  =  ic/2  and  p  gates  from  0  to  -00,^6  quantity. 
njS-sin  a  slnh  (3  Is  Initially  a  decreasing  function  of  p. 

The  integrand  above  has  singularities  at  p^  =  p  cos0c  + 
ne  +  n(2kn)  if  sine  =—  (k  =  0,  +  1,  These  are 

O  Op  ““ 

arranged  on  the  p-plane  in  a  sequence  parallel  to  the 
Pj  -  axis  at  intervals  of  k(2nn);  we  will  show  in  a  moment 
that  the  are  branch  points.  Thus  the  contour  X  encloses 
only  a  single  branch  point;  it  is  this  circumstance  which 
determined  the  choice  of  C* ,  for,  asusual,  we  will  wish 
to  apply  Cauchy’s  theorem  and  close  about  a  vertical  branch 
cut , 


Consider  p  =  pcos0  +  n0  ;  we  want  the  behavior  of 
01  o  o 

n  -  p  sine  in  a  neighborhood  of  pQ,  that  is,  when  0  approaches 
0Q.  Upon  expanding, 

n  -  P  sin©  =  -  p  cos0 .'(0-0  )  +  (  Psin0o)(e-e  )2  +  ... 

I  I  e  o  —  o 

p  =  P  cos©  +  ne  +  (-  P  sine  +  n)(e-e  )  +  ~ 

|  o  o  |  0  o 


*  po  “  ^-C-gy9°-  (0~eo)2  +  ••• 

Therefore,  near  p.,  (p-p_)  behaves  like  -  Pcos0o(e-e  )2. 

•  r0  O  — g  '  »■  *  o  9 
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and  so 

n  -  p sin©  /-v/  - 


-,1/2 

2(P“PQ}  |  +  0(p-p„)  + 


cose. 


This  indicates  that  pQ  (  and  p^  for  all  k)  are  branch 
points.  The  contour  contains  pQ;  we  take  the  branch  cut 
vertically,  and  find  that  after  shifting  to  Y  1 , 


>  I  •  f 

f  '1 

I  1 


|  r-;( 


.  V 


,oo 


i(p0+is) 
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(  ^cose  ) 


ids 


Here  p  =  pQ  +  is  and  the  first  term  in  the  expansion  of 
n  -psine  has  been  osed  for  the  asymptotic  estimate*  Thus, 

Hn1}(p  \jfni/2)  eIp° . 

Now,  £  =  sineQ  *  sina0cosh(3o  +  I  cosaQcoshPo  -*oo  as  n  be¬ 
comes  large, Because  jr  is  real  foyr  .any  value  of  n,we  have  aQ=it/2: 

referring  to  the  contour  C*  in  the  (a  +  ip)  =  ©-plane,  evi- 

p  -P  -P 

dent  ly  aQ  =  n/2  ■ — y'  p^  <  0,  Thus  j-  —  ^(e  0  +  e  *ye  ■  ? 

for  n  lar&e,  implying  that  PQ  -ln(^~)  , 

Also, 

pcoseQ  *  ^U-sin2©o)1/2  =  (  p2  -  n2)  */2  , 

p0  =  pcoseo  +  n©o  =  (f2  -  n2)1/2  +n(a0  +  ipQ), 
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which  indicates  that  as  n  becomes  large, 

1 3  cos0o  ^ 

P0  r-s/  in  -  in  in  (^p)  +  . 

Consequent ly, 

„U)(r  ,  ^  IdL Jlfj  rufi,  JTn  e1Un  +  ‘  lnIn‘^ 

or,  _ 

)'v  j|  <nntl/2e -")(£)*"  . 
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XI.  Field  produced  by  a  source  in  the  presence  of 
a  conducting  wedge* 

For  simplicity,  the  two  dimensional  situation  will  be 
considered.  Assume  the  wedge  to  be  of  infinite  length  with 
its  edge  along  the  z-axis,  and  let  a  sinusoidal  line  source 
parallel  to  the  z-axis  be  given  at  (rQ,  0o). 


The  electromagnetic  field  components  will  be  independent 
of  z.From  Maxwell's  equations  (eq.26)  T.R.  No. 2  it  is  easily  sc 
that  two  possibilties  exist:  we  may  either  solve  for  Ez  from 
which  will  follow  H  and  H  ;  or  solve  for  H  which  also 

a  y  6 

gives  E  and  E  ,  Both  E_  and  H  are  seen  to  satisfy  the 
a  y  z  z 

wave  equation 

=  ■  Mr-rn)  £(9  -  eo) 

z  r 

Supposing  the  wedge  to  be  a  perfect  conductor,  the  boundary 
conditions  are; 

( i )  Ez  =0  at  9=0  and  0  =  2k  -  (3 

( i i )  “  Hz  =  0  at  0  =  0  and  0  =  2k  -  0  . 

■De 
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Although  two  distinct  problems  have  been  formulated 
the  methods  of  solution  are  essentially  identical.  We  shall 
solve  for  E  and  in  the  course  of  the  argument  indicate  the 
differences.  Of  particular  interest  in  the  solution  *  • 
will  be  the  shadow  region  it  +  0Q  <  9  <  2%  -  (3, 

whose  “darkness"  varies  as  a  consequence  of  diffraction  at 
the  apex  of  the  wedge.  Two  other  interesting  situations 
are  the  limiting  cases  (3  *  %  and  p  =  0;  the  former  has 

already  been  dealt  with  in  Section  IX  (reflection  off  a  plane 
conducting  surface),  while  the  latter  is  the  condition  for 
the  Somme rf eld  half -p la ne  problem  (diffraction  by  a  knife  edge) 
The  Sommerfeld  problem  will  be  discussed  below. 

As  in  the  preceding  section,  the  wave  equation  may  be 
written  in  plane  polar  form: 


I  3L  + 

r  ^  r 


+ 


j(f-r0)^(e-e0) 


The  equation  is  separable;  therefore  put  ••  *  L,  L  is 

2e2 

an  operator  acting  on  the  space  of  functions  u(9), 

0  <  9  <  2n  -  p,  which  satisfy  the  boundary  conditions  (i) 
above, The  eigenfunctions  of  L  are  derived  f ron(  ^  I  -  L)u  *  0; 
they  are  sin  ^  0,  with  ^  •  (Were  we  working  with 

H£  the  only  difference  here  would  be  that  u(e)  satisfies 
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the  boundary  conditions  (ii),  necessitating  the  use  of 
cosines  instead  of  sines,)  We  have 


,  d1  ,  i  a  .  k2  l,  .  kr-r0U(®-V 

(i7*7-S?  +  lt  '7)e* - r 

The  homogeneous  equation  is  Bessel's  equation  (eq,  2;');  as 
in  Section  I,  solutions  must  be  chosen  which  are  regular 
at  the  origin  and  outgoing  at  infinity.  These  were  labeled 
J^(Kr)  and  H^fkr),  respectively,  so  that 


Ez =  J\rc{kr)H^)(kro)a  r<r0 

*  JyfI(kro)H}~)(kr)a  r  >  r0 

where  a  is  the  Jump  factor,  Reference  to  Section  1  gives 

a  =  -  j  however,  note  that  in  Section  I  we 

found  the  jump  f  actor  by  applying  the  asymptotic  forms  of 

the  Bessel  and  Hankel  functions  to  the  Wronskian,  and  that 

these  forms  were  derived  for  integral  orders  of  J  and 

r  r 

only.  It  will  be  shown  in  a  moment  that  all  the  previous 

representations  of  Jr  and  are  valid  for  non-integral  r, 

which  will  Justify  a  *  -  J(9-0O)  • 

To  interpret  the  order  /L  of  the  Bessel  and  Hankel 

functions,  we  expand  <$(0-0o)  in  eigenfunctions  of  L,  Letting 

so  that  ^  =  n  Y  #  it  can  be  shown  that 
co 

(T(0-9o)  *  ^  sin  ny0  8innVeo  , 
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2  v  # 

where  —■  Is  the  normalization  factor.  The  operator  L  appearinj 
It 

in  E„  can  now,  as  in  Section  VIII,  'Tech,.  Report  No.  2,  be  ‘ 

2  ,2 

replaced  by  its  eigenvalues  nj"  .  Thus 
(  1$)  Ez  =  1  X  &  Jn  kr<: )H^ ^ ( kr> )  sinn(f ©sin  n  2f  Qq  ; 

r>  -  max(r,rQ)  and  r<  =  min(r,rQ)  has  been  used-.  (H^  is 

given  by  equation  (l5)upon  the  replacement  of  sinn^®  by 

wo 

cos  n)f  ®  ,) 

0o 

In  order  to  check  the  convergence  of  Ez  we  utilize 
approximations  of  and  which  are  valid  for  large 

values  of  A.  With  recourse  to  the  familiar  series  expansions, 

JAP  )~  j  H(1)(o)  ^P(A)(  P/2)’A,  A  >  >  1* 

A  ’  r  (  +D  ^  r 

That  is,  when  A  is  large  and  behave  like  the  first 

terms  in  their  expansions.  Therefore,  for  sufficiently  large 

n» 

Jn  *  (kr<,Hi1,)(kr>)  ~  • 

Since  ^  <  1  wherever  r  f  rQ,  we  conclude  that  the  series 
in  equation  (7)  converges  absolutely  except  on  the  cylinder 
r  =  r0,  (Although  of  no  physical  significance,  it  is  possible 
to  prove  conditional  convergence  at  r  =  rQ,  a  fact  due  to  the 
presence  of  the  sine  functions,) 

The  behavior  of  J N(  p)  at  small  values  of  p  is  again 
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dominated  by  the  initial  term  of  its  series  expansion. 

Consequently,  near  the  apex  of  the  wedge  we  may  approximate 
oo  ny 

Ez  ^J  l*  2  Hn1)(kr0)8in  **  6  n/  9Q  ^ 

1 

'V'  H^(kr0 )  sin^ e  sin/  0Q,  kr  <  <  1. 

This  approximation  indicates  the  qualitative  behavior  of 

E  (  and  H  )  at  the  wedge  apex.  Consider  the  limiting 
z  z 

cases  p  «  o  and  p  =  ®,  i.e«,  ^  *  1/2  and  =  1: 

For  P  =  0,  E  =  0(r*/2)  and  H  and  H  ,  which  are  propor- 
z  x  y 

tional  to  the  derivatives  of  E  .  are  OCr-1^).  Thi  s  shows 

z 

that  the  Poynting  vector  E  x  H  is  bounded  as  r-0  and 
therefore  the  energy  of  the  field  is  finite  in  the  vicinity 
of  the  knife  edge  even  though  the  magnetic  field  is  infinite 
there. 

For  p  =  it,  E  *  0(r)  and  then  H  and  H  will  be  0(1). 
z  x  y 

This  shows  that  both  the  electric  and  magnetic  field  are 
finite  near  r  =  0.  This  result  was  to  be  expected  since 
the  wedge  in  this  case  reduces  to  a  plane. 

We  now  turn  to  an  Investigation  of  the  field's  beha¬ 
vior  in  general.  This  will  become  more  feasible  if  the 

infinite  sum  representing  E_  can  be  replaced  by  a  more 

z 

manageable  expression;  and  as  a  matter  of  fact,  upon 
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+  Ik2  )  )  eikr  cos  <<P-«/S>dq,  . 

7  C1 , 

( 

Consider  -  -j  J  s 

r  q 

^  J  .  eikr  cos  <pei/((p-n/2)| 


7  ^  "  ‘7 

r  C,  r 


y  tci 


r 

ikr  sin  <p  eikr  cos  ^e1  ^<P"V2) 


dtp  = 


1  ikr  sin  ©  _ikr  cos 
rr/ - >r- ,  — *  e 


7  t?-- 


+  ... 


...7$  5 


2  2  2 

k  r  sin  tp  -  ikr  cos  cp)e 


ikr  cos<pci^(<p-ii/2 

~Ti 


when 


V  2  \ 

- — y  is  substituted  into  [  ] J  ,  we  find 


[  ]  3  =  integrated  part  « 

ci  L 


-  / 1 _ ikr  sin  P\  ikr  cos  <p  if  (<p-it/2)l 

Y2  /  e  lci  • 


It  is  not  difficult  to  show  that  the  integrated  part  is 
equal  to  zero*  The  integral  representation  of  is 

therefore  valid  for  all  real 

The  computation  of  the  asymptotic  form  of  followed 

from  the  integral  representation;  the  argument  of  pages  12-16 
with  minor  modifications,  will  also  go  through  for  non- 
integral  $  .  Referring  to  equations  1$  , 
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HyWw  ^je1{kr  ”  ^  Kr  >  >y. 


Since  H^1}  =  h(?^,  we  find  as  before 

* 


H^J(kr 


)  =  1  ^  eIkr  cos  ‘Pe  (‘P'lt/2)d<p. 


c2  (- 


*}.  _  ni  (xkjl. 


T  5T r 


also. 


0  <  <  cu,  <  2n,  joj-agl  =  ft 


H(2)(kr)^->  vLjpe~i(kr  “  “2f  */2)»  kr  >  >  ^ 


The  Bessel  function  Jy  may  still  be  expressed  as  the 

fl 

following  linear  combination  of  the  Hankel  functions: 

(kr)  =  2  Hy  }(k«F)  +  £  H(2,}(kr)  . 

But  is  required  to  be  regular  at  the  origin  and  this 
must  be  verified  for  the  above  sum  when  is  not  an  Inte¬ 
ger,  Upon  use  of  the  integral  representation,  it  is 


evident  that 


iV  (kr)  =  -n \rS 


(Ikr  cos  <pe  i2T(q>  -  */2  )d(j)  * 
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Thus, 

yo) 


elf  (<?-V2)d<p 


which  for  >  0  is  easily  seen  to  be  finite. 

Although  the  integral  expressions  for  J  and  are 

now  available  for  use  in  equation  15  ,  E  at  best  remains 

an  infinite  sum  of,  rather  than  a  single,  contour  integral. 

To  write  Ez  in  the  desired  form,  and  thereby  proceed  with 
the  examination  of  its  general  behavior,  the  original 
situation  will  be  slightly  simplified.  Instead  of  consi¬ 
dering  a  cylindrical  wave  incident  on  the  wedge,  it  will 
be  assumed  that  the  radiation  source  is  a  plane  wave  in¬ 
coming  from  the  direction  0Q.  In  effect,  the  assumption 
is  tantamount  to  letting  the  original  line  source  go  to 
infinity  (rQ  -*  oo).  Though  the  source  is  altered,  diffraction 
still  occurs  at  the  wedge  apex,  the  region  k  +  Qq  <  0  <  2n-{5 
continues  to  lie  in  the  shadow,  and  a  wave  is  still  reflected 
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off  the  surface  of  the  wedge.  The  incident  plane  wave 

case  therefore  presents  a  good  —  and  solvable  —  analogue 

to  the  original  wedge  problem. 

We  know  the  series  of  equation  l£  is  absolutely 

convergent  for  r  4  r  .  Since  r  -*•  oo  .  E  may  be  written  as 
'O  0  r  z 

N  _ 

Ez  =  llZ:  (*T>  fife  n/esl„ 

from  which  it  will  follow  that 

Jn«<kr)'IS?  ,“kr0-«A-»ir  V6),ln  „jre  5l„  ntfe0  - 
1  0 

a  e1(kr°-U/^)^eikr  cos  *  ^  ein^(<p‘ll)sin  n)fo  sin  nJeoc 


As  r  -*  go  ,  |  El  «•  0;  that  is,  the  amplitude  of  the  field 
falls  off  to  zero  as  the  source  is  moved  to  infinity.  Thus, 
also,  the  incoming  part  of  E_,  which  will  be  the  incident 
plane  wave,  loses  its  amplitude  as  r  oo.  In  order  to 

r«kr"fi  “i(kr 0~*A) 

nil  -  1  A  f  k  1  p  •  1  A  1^.1--  r  (  I  jtA*  ^  V  #  _ 


avoid  this,  we  multiply  E  by 

2 


—  in 


effect,  E,  is  normalized.  The  z-component  of  the  electric 
field  in  the  presence  of  the  wedge  is  therefore 


F  =  \ 


fitkr"  -i(kr  -n/U 
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We  get 


F  elkr  cos  t  Z  ein^(<P“*)sin  „j9  ,ln  nfe^ 


00 


because  we  may  write 


sin 


w  ^  ,  fin^e-eJ  -in/(e-ej  in2f(©*e)  -in2f(©+0  )1 

n^esdn  ntfeQ=  £je  °+e  0  -e  0  -e  °J, 


F  becomes 

FJi  (  ikr  cos  cp  >  C(1^((p”lt+e“e o))n.  e 

F"Hi  l  e  4>Le  +*““e 

3 

v  (  i  i  ((P-H+0-0  ) 

-  ki  )  e*kr  »s  »,  e  '*  o'  +___. 

otc  ■(.  - ! — i  _  — - „■  -r  ••• 


( i  ^(cp-it-e-e  ))n 


n  I 

Jdcp  , 


i  ( cp-'Jt-e-e  ) 


i  (<p-ji+9-0 _ y  i  (cp-it+6-9  ) 

1-e  °  1-e  ° 


dtp 


Thus  the  field  incident  on  the  wedge  has  been  represented  by 

a  single  contour  integral.  The  sums  of  the  first  two  terms 

in  the  square  brackets,  and  the  second  two  are,  respectively 

cos^e-e  )  -  e*  <p— it )  cos/(0+0  )  -  e  *  ^ ) 

CO  Sjf  ( cp-u )  -CO  s£  (  0-0Q )  "  cos^ltp-it)  -  cosn0+0o)  • 

When  <p-K  is  replaced  by  cp  and  the  contour  changed  accordingly, 

we  have  after  addition, 

slwfo 


F  =  £)  e“ikr  C0S  * 

°k 


Hi 


I, 


COS  <p-COS  (e 


=*] 


dtp 


ikr  cos  cp 


cm  ;a-u0l]  **  * 


F1  *  F2* 


-  3*  - 


rs  t. 


-IT  -«T 


-2n  <  <  a2  <  %/  (a^-Ogl  = 


Although  neither  F^  nor  F^  can  be  explicitly  evaluated, 
we  can  nevertheless  examine,  as  usual,  their  asymptotic  be¬ 
havior  and  thereby  arrive  at  significant  results.  Consider 
Fj  for  ((3  *  0  ^  -  1/2 j  without  significant  loss 

of  generality,  we  may  reduce  the  wedge  to  a  knife  edge).  Then 

(  sin  2<9 

p  _  _ 1  )  -ikr  cos  <p  ■■  ■  — ■!'"■■■  -  ■  ■  ■  i ■  —  j 

1  Th%  g  COS  £(p  -  cos  ^  ( 0-0Q )  ^  * 

Let  ^  =  cos  To  ~  cos  ^(9“0O )»  d?/ »--iysin  ^q>  d<p.  Notice 

that  the  zeros  of  sin  ^  <p  determine  branch  points  in  the 
X  -plane,  i.e.,  cos  j  2nx  **  +  1  are  branch  points.  Thus, 


p  _  2  \  -ikr._. _2 

3  t  -  x. 


*T  =  -^re^r  j  e-***r  y* 
B  r-r, 


<*r , 


where  B  is  to  be  computed.  The  contour  may  be  conveniently 
taken  with  *  -  3it/2  and  a2  =  */2  . 

a)  When  <p  “  -  +  ip,  ^  **  Tj  +  i?^  *  -cos  ^  cosh  Jp  + 

+  i  sin  ^  sinh  ^  *  "tan  taifc  P/2  "*  "l 


as  p  -*  00  j 


the  contour  B  •*.  approaches  a  slope  of  -  1 


In  the  second  quadrant.  p  *  0  ^  -  1/  /2~. 

b)  When  <p  *  a+l£  (  thesis  to  prevent  B  from  passing 
through  the  pole  ^ )  : 

-  3it/2  <  a  <  -it  ^  i  <  0  , 

“*  ”  ¥2  1  “* 

-  n  <  a  <  0  ^  0  <  ^  <  1  , 

0  <  a  <  n/2  1  >  7^  >  1//T  , 

and  may  be  put  equal  to  fc'on  this  segment  of  B. 

c)  When  <p  **  n/2  +  ip,  -  -1  as  p  -*  oo  ,  so  that  B 

X  1 

also  approaches  a  slope  of  -  1  in  the  fourth  quadrant. 


We  wish  to  shift  the  contour  B  so  that  e“2ikr£ 
becomes  an  exponentially  decreasing  function.  The  line  1  . 
in  the  diagram  is  suggestive;  and  in  fact  on  this  line 
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*2/  s»  -  oo <  R  <  00  ,  and  *  -iR2,  which  is  just 

what* s  required.  The  loop  about  the  branch  point,  however, 
presents  a  -problem,  Cauchy's  theorem  cannot  be  applied  to 
the  region  containing  +1,  and  so  B,  after  it  is  shifted  to^  , 
will  have  to  be  left  with  a  loop  about  +1,  This  difficulty 
can  be  neatly  avoided  by  noticing  that  the  integrand  of 
is  an  odd  function  of  <p,  which  means  that  the  contributions 
from  -  it/2  <  9  <  0  and  0  <  <j>  <  it/2  cancel  each  other.  It  is 
just  this  Interval,  -  n/2  <  (p  <  h/2|  which  has  given  rise 
to  the  loop  in  the  f -plane.  Hence,  by  neglecting  it,  it 


The  pole  will  have  to  be  taken  into  account,  that  is 

when  B  is  shifted  to  £  we  will  acquire  a  residue  if  V  >0, 

**  o 

The  residue  at  t0  is 


■2ikr  X2 
wo 


2 

i  -ikr(2^0  -1) 

r 


*  i«"Ikr  cos(e-e0) 


The  change  of  sign  is  due  to  the  orientation  of  B  in  the  lower 


-  je  - 


half-plane.  Thus, 


ilt,.  \  o 

F1  “  dh  +  H(  V  )i  e-ikr  «•  (6-80)  > 

d  y  -  c0  0  ^  ' 


where  H(^q)  =  1  when  >  o,  aero  otherwise.  Or, 


when 


f2  =  2krr2,  .A 


-eikr  f  ,*2  1 

Fi  =  tt  \  *  J  jZrM? 6i  +  H(*o>!Klkr  cos(e-eo> . 

■rl  -iti/k  0 


In  order  to  estimate  the  integral  in  F,,  we  use  a 
technique  made  familiar  in  preceding  sections.  Set 


3-^/STr  \J(Sr  X0  '>2kr  '  i  -  y  fm  ; 


then. 


•ni/k 


>~*lA 


-jCl  \  e"1*  d  6  cIkr 


\  e-^2_ 

-BIT  \  i 


**-‘3 


—u>  e 


-iti/V 


'*  )7dl*RT  Vj-’e'o&trl 

4o  e^A 


el(kr-TC/lj.) 
8  ^  \j2itkr 


♦  Ei  » 


where 


-  3‘9  - 


>Eii  *  -fe  )  e'R2 1  R^«iAr0^  I*  •  r.isr  • 


The  minimum  value  of  the  denominator  in  the  integrand  of  the 
error  term  is  simply  the  distance  of  e71*^  \)  2kr  from  the 


real  axis,  namely 


1/  \fkr". 


Theref  ore 


oc 

|£>'  £  8 j 


|R|e-R  dR  =  0(— • 
krt0 


Ej  is  a  satisfactory  error  term  so  long  as  ^  is  not  in  the 


neighborhood  of  zero*  We  have 


.  -ikr  cos ( 0-0  ) 
^-HlVo^e  0  - 


i(kr-xA) 


QT?  42n kr 
o 


Ej  , 


Similarly  for  F2;  if  'fc'0  =  cos  £(9+0o)  * 

F,  rn  JI(VI  )ie"lkr  c°s(e+0o  )  -  e1^kr 

2  u°a  2 

0 

The  field  produced  by  the  incoming  plane  wave  in  the 
presence  of  the  knife  edge  is  Fj  +  F2  =  F: 

_  .  -ikr  cos(e-e  )  ,  -ikr  cos  (0+0  ) 

(16)?=  H(to)^e  0  -  H(2T^)|e  0  +  ... 


...  + 


i(kr-*/U)  /  v 

d\l2vkr  \  'to  *t>  q» 


,)  +  E- 


8 1Z5r  Wo  ' 

where  7/1  ■  cos  £  (0-0Q),  *£'  *  cos  £  (0+0  )  » 


ko 


(^+  *?)]• 


The  physical  interpretation  of  equation  (16) is  clear, 

.  -ikr  cos(e-0Q) 

The  first  two  terms  of  F  are  plane  waves;  ij^e  u 

is  the  source  wave-  incoming. along  the  direction  0Q,  and 


-ikr  cos(e+eJ 


is  incoming  along  -0Q,  i.e,,  outgoing 


along  it-0Q,  and  it  represents  the  reflection  of  the  source  by 

the  knife  edge.  Since  V  >  0  when  0  <  0  <  it  +  0  and 

o  o 

V 

u  *Q  >  0  when  0  <  0  <  it-  0Q,  we  know  in  which  regions  these 
plane  waves  exist.  The  remaining  term  in  equation  (l6)is  a 
cylindrical  wave  due  to  diffraction,  and  it  exists  everywhere, 

.IT-*.  /ft 


m,  :  £  *o 


Jl 


In  I,  F  *  inc.  +  outg,  +  cyl,  *  source  +  reflec.  +  diff,* 
in  II,  F  =  inc.  +  cyl,  *  source  +  diff,  ; 
in  III  (shadow  region),  F  «  cyl.  =  diff. 

Equation  (16)  yields  no  information  about  the  behavior 


of  the  field  along  e  *  x-  0„  and  0  *  n+  0^  since  there  either 

oo 

or  equals  sero  and  the  error  term  breaks  down. 
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In  order  to  determine  the  field*#  behavior  on  the  boundary 
lines,  we  will  have  to  arrive  at  an  estimate  of  F  which 
does  not  depend  on  the  magnitude  of  V  or  V  * 

O  <^Q  * 

Consider  Fjj  from  before, 

CO  .-"‘A 

COSle'9o)  •  &  ( 

-<o  e”*‘A 

tet  1  -  to-«‘A  Md  Rq  -  e"Av  .  th.  Integral 


Fj  becomes 


I  .  -o^r 

1  ~sr 


oo 

(  .** 

)  ttr**  • 


Since  R-Rq  *  R  -  X  QVkr  -  i  U  v/kT  , 


we  can  write 


.1  (  Ir(R-R^) 

TOnT)  *  \  •  °  dr  if  XQ  <  0, 


i  (  -ir(R-R^) 

TTOT)  “  \  e  dr  if  X  >  0. 

O  )  o 


Therefore,  supposing^  <  0,  we  have 


F1  *  ll  *  l£S«"  J  «“R2dR  |  e 


Ir(R-Ro) 


-  ki  - 


A  similar  result  obtains  for  Fg  :  in  the  neighborhood 
of  0  =  ft-0Q,  ^2  ^  H  Notice ‘that  both  and  are 
outgoing  along  the  regional  boundaries. 

With  reference  to  eqaation(i6)  then,  we  have: 


near  0  =  n-0o  ^  0):  F  ^elkr  +  £  e“ 

i(kr  -  n/k) 


ikr  cos (0-0  ) 


V  8j2*kr 


fe) 


near  0  *  *+0q  (£ Q^>  0):  F  ^  je 


i  ikr  e 


i(kr-nA) 


X  '8V5i3tr 


>(-4) 


In  these  expressions,  ^  ^  -  cos  0Q  and  VQ  ^  cos  0Q, 

Equation  (16)  gives  the  field  for  all  other  values  of  0, 


